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Abstract 



The main general relativistic effects in the motion of the Moon are briefly 
reviewed. The possibility of detection of the solar gravitomagnetic contributions 
to the mean motions of the lunar node and perigee is discussed. 
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1 Introduction 



In a recent paper, Gutzwiller has provided an admirable review of the oldest 
three-body problem, namely, the Sun-Earth-Moon system [1]. Some work on 
the relativistic theory is mentioned in his paper; however, in view of the recent 
advances in relativistic celestial mechanics this subject deserves a more complete 
discussion. Here we provide a brief description of the main relativistic effects. 

The lunar laser ranging experiment has opened up the possibility of mea- 
suring relativistic effects in the motion of the Moon; indeed, the agreement 
between the standard general relativistic model that contains over a hundred 
model parameters and the ranging data accumulated over the past three decades 
is excellent [2,3]. For instance, the post-fit residuals in the Earth-Moon distance 
are at the centimeter level [2,3]. Simple theoretical estimates lead to the con- 
clusion that the main relativistic effects in the lunar theory are due to the 
spin-orbit coupling of the Earth-Moon system in the gravitational field of the 
Sun. The post-Newtonian influence of the solar field on the lunar motion con- 
sists of terms that can be classified as either harmonic (i.e. periodic) or secular 
(i.e. cumulative) in time. It turns out to be very difficult in practice to separate 
the harmonic terms from the corresponding Newtonian terms with the same pe- 
riodicities. In effect, the existence of the post-Newtonian harmonic terms leads 
to small relativistic corrections in the numerical values of certain model param- 
eters that are thereby adjusted by a fit to the ranging data. To give an example 
of such harmonic effects, we mention our prediction of a 6 cm relativistic tidal 
variation in the Earth-Moon distance with a period of 1/2 synodic month [4]. 

The main secular terms turn out to be essentially due to the precessional 
motion of the Earth-Moon orbital angular momentum in the field of the Sun. 
The Earth-Moon system can be thought of as an extended gyroscope in orbit 
about the Sun; we are interested in the description of the motion of the spin 
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axis of this gyroscope with respect to the "fixed" stars (i.e. the sidereal frame). 
An ideal pointlike test gyroscope carried along a geodesic orbit would exhibit, 
in the post- Newtonian approximation, geodetic precession due to the orbital 
motion around the mass of the source as well as gravitomagnetic precession due 
to the intrinsic rotation of the source; however, the finite size of the gyroscope 
in this case (i.e. the orbit of the Moon about the Earth) leads to additional tidal 
effects. The main post-Newtonian gravitoelectric effect, i.e. geodetic precession, 
results in the advance of the Moon's node and perigee by about 2 arcseconds 
per century as first predicted by de Sitter already in 1916 [5]. This motion has 
been measured by Shapiro et al. with an accuracy of about one percent [6]. 
It is a relativistic three-body effect; therefore, we consider in the next section 
the restricted three-body problem in general relativity and briefly indicate, in 
particular, the more subtle post-Newtonian gravitomagnetic contributions to 
the motions of the Moon's node and perigee that are caused by the rotation of 
the Sun; indeed, solar rotation induces cumulative relativistic tidal effects in the 
Earth-Moon system [4] . 

2 Restricted Three-Body Problem in General 
Relativity 

In our previous work [4] , we developed a new scheme for the approximate treat- 
ment of the restricted three-body problem in general relativity. This coordinate- 
invariant approach is particularly useful for a reliable theoretical description of 
relativistic (solar) tidal effects in the motion of the Earth-Moon system. We 
assume that the Moon follows a geodesic in the gravitational field generated by 
the Earth and the Sun. This field may be calculated as follows: we first imagine 
that the Earth follows a geodesic in the solar field. Along this geodesic, we 
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set up a geocentric Fermi coordinate system. This system, which involves the 
tidal field of the Sun, is then enhanced by taking due account of the field of 
the Earth in the linear approximation. Tidal effects in general relativity involve 
the projection of the Riemann tensor onto the tetrad frame of the measuring 
device. Consider the tidal matrix for a test system ("Earth") in free fall in 
the gravitational field of a rotating mass ( "Sun" ) . In the standard first-order 
post-Newtonian treatment, the spatial axes of the local tetrad frame along the 
orbit are obtained by boosting the background Minkowski axes and adjusting 
scales to maintain orthonormality; the resulting tidal matrix for an approxi- 
mately circular geodesic orbit turns out to be sinusoidal in time [7]. In this 
case, the tetrad frame is not parallel-transported, but its motion involves the 
Lense-Thirring orbital precession as well as the geodetic (i.e. de Sitter-Fokker) 
precession of the spatial axes. Once the parallel transport of the spatial axes 
along the orbit is imposed, the gravitomagnetic (i.e. Schiff) precession of the 
spatial axes would also appear in the first post-Newtonian order. In this or- 
der, the tidal matrix for the parallel-transported axes contains a secular term 
as well that must therefore be a direct consequence of the Schiff precession of 
the spatial axes [8], in agreement with our previous work [9-1 1]. The linear 
growth of this gravitomagnetic contribution to the tidal field poses a problem 
for the first post-Newtonian approximation: the non-Newtonian "off-diagonal" 
part of the tidal matrix can diverge in time [9-11]. To avoid this limitation, 
we have developed a post-Schwarzschild treatment of gravitomagnetic tidal ef- 
fects; indeed, the concept of relativistic nutation provides a natural resolution of 
this difficulty by limiting the temporal extent of validity of the post-Newtonian 
approximation [9-11]. 

Imagine, for instance, a set of three orthogonal test gyroscopes falling freely 
along an inclined circular geodesic orbit with constant radius r ("astronomical 
unit") about a slowly rotating central body ("Sun") with mass M and proper 
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angular momentum J. The motion of the spin axes of these torque-free gyro- 
scopes, which constitute a local inertial frame (i.e. the geocentric Fermi frame), 
is essentially governed by the equations of parallel transport along the geodesic 
orbit. By solving these equations using the post-Schwarzschild approximation 
scheme that takes M into account to all orders, it can be shown that the av- 
erage motion of the gyroscope axes with respect to an effective Newtonian (i.e. 
sidereal) frame consists of a gravitoelectric precessional motion — i.e. geodetic 
precession that was first completely analyzed by Fokker — together with a com- 
plex gravitomagnetic motion that can be loosely described as a combination of 
precessional movement and a harmonic nodding movement. The latter motion 
is a new relativistic effect of a rotating mass and has been referred to as relativis- 
tic nutation [11]. In the post-Newtonian approximation, the nutational terms 
over a limited time combine with the other gravitomagnetic precessional terms 
to give the Schiff precession. To see how this comes about, let us denote by r 
the proper time of the geodesic orbit and consider a vector normal to the orbital 
plane (ecliptic) at the beginning of measurement (r = 0). Rclativisticnutation is 
a periodic variation of the angle between this vector and a gyroscope axis that 
is Fermi propagated along the orbit. The leading contribution of relativistic 
nutation to this angle can be written as 

9 n w £[sin(?7o + ujft) - sinryo] sin a, (1) 

where T]o is the azimuthal position of the Earth in the ecliptic at r = measured 
from the line of the ascending nodes and £ = J/Mr 2 uo. Here uj,oj 2 = GM/r 3 , 
approximately describes the orbital frequency in the absence of rotation and a 
denotes the inclination of the orbit with respect to the equatorial plane of the 
Sun [12]. The frequency of this nutational oscillation is the Fokker frequency 
oj f rs few, where e = GM/c 2 r. The nutation amplitude, £sina, does not 
depend on the speed of light c. This remarkable fact can be traced back to 
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the occurrence of a small divisor [9-11] involving the Fokker frequency. In the 
post-Newtonian limit of the post-Schwarzschild approximation, Eq. (1) reduces 
to 6„ <~ u) n T, which represents a precessional motion with frequency oj n = 
£u>f sinacos^o about a direction opposite to that of orbital velocity at r = 0. 
Thus, relativistic nutation reduces to a part of the Schiff precession in the first 
post-Newtonian approximation. It follows from this analysis that the first post- 
Newtonian approximation breaks down over timescales of the order of Fokker 
period tf — 2it/ujf; however, this fact does not diminish the usefulness of the 
first post-Newtonian approximation for the description of observations in the 
solar system since in this case the Fokker period is almost immeasurably long 
(e.g. tf — 67 million years for the motion of the Earth about the Sun). 

Let us consider the influence of the gravitomagnetic field of the central body 
("Sun") on the relative acceleration of two nearby test particles ("Earth" and 
"Moon") moving along the circular geodesic orbit. The dominant contribu- 
tions of the gravitomagnetic field of the central body to the tidal matrix, first 
calculated by the authors [9-11], are proportional to 



which is directly proportional to the amplitude of relativistic nutation (£ sin a) 
and exhibits a maximum (at t = tf/2) that is independent of the speed of light 
c. It follows from Eq. (2) that to first order in loft <C 1, the dominant gravit- 
omagnetic amplitude varies linearly with r. This secular amplitude originates 
from a coupling of the nutation part of Schiff precession with the amplitude 
(~ uj 2 ) of the Newtonian contribution to the gravity gradient [13]. It should 
be mentioned in passing that the relativistic quadrupolc contributions to the 
tidal matrix have properties quite similar to the gravitomagnetic tidal effect 
described here [10]. 

Let us now turn to the potentially observable effects of the solar gravit- 
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omagnetic field on the lunar motion. The lunar path is determined by the 
Newton- Jacobi equation 

d 2 x i Gm ; ■ 

- f - T + 1 pX* = -K ij {T)x>, (3) 

where x l ,i = 1,2,3, represent the geocentric Fermi coordinates of the Moon, 
m is the total mass of the Earth-Moon system and R(t) denotes the Earth- 
Moon distance depending on the proper time r measured along the geocentric 
path around the Sun. Here K is the tidal matrix. Equation (3) describes the 
motion of the Moon with respect to a geocentric local inertial frame [14]. Using 
the equation of relative motion (3), we have calculated — among other things — 
the influence of the tidal field of the Sun on the orbital angular momentum 
of the Moon with respect to the Earth. To express the result with respect 
to the sidereal frame, we choose as our sidereal reference framethe geocentric 
Fermi frame at r = 0. This frame is related to the Fermi frame at time r by 
a rotation matrix that incorporates the relativistic precession and nutation of 
the Fermi frame with respect to the sidereal frame. In the first post-Newtonian 
approximation, this motion reduces to a (Fokkcr plus Schiff ) precession. Let D 
denote this rotation matrix, then d = DijLj, where the sidereal components 
(£i) of the orbital angular momentum are obtained from a transformation of 
the geocentric components (Lj) with 

Dij = Sij - e ijk <f> k , * = / uj FS {T')dT'. (4) 

Jo 

Here the analysis is limited to the first post-Newtonian approximation and ujfs 
represents the frequency of (Fokker plus Schiff) precession. The direction of 
Schiff precession is not fixed along the Earth's orbit; therefore, $ contains (cu- 
mulative) secular terms (which represent simple precession) together with (har- 
monic) nutational terms of frequency 2ui and amplitude of order aet;. Averaging 
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over the latter terms, the dominant secular terms in <fr are given by 

$1 ~ -ae£ojT sinr^o , &2 ~ ^ — 2 ^) tL ° T ' ^ 3 ~ 2 Q;e ^ a;TCOSr ' ^ 

with respect to the geocentric Fermi frame at r = 0, which has its first axis 
essentially along the radial position of the Earth, the third axis approximately 
along the direction of motion of the Earth and the second axis normal to the 
ecliptic (in a direction opposite to the Earth's orbital angular momentum about 
the Sun). We note that for £ = 0, Eq. (5) expresses the de Sitter- Fokker effect 
that has been observed by Shapiro et al. [6]. To illustrate our approach, let us 
use Eq. (3) to determine the value of L, which is the angular momentum of 
the Moon in a circular orbit about the Earth with respect to the Fermi frame, 
averaged over orbital motions of the Earth about the Sun (with frequency oj) 
and the Moon about the Earth (with frequency SI). Then 

f^^xLo , (6) 

where Lo is the unperturbed orbital angular momentum with respect to the 
geocentric Fermi frame and Cj is given by 

u) 1 ~ — tl> ae£ ^2 sin 770 + ^wt cos 770^ , (7) 

^ 2 «^ (l-6e£), (8) 



OJ 



wo«£( ( 2cos?7o — —uJTsinrjo ) (9) 
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to first order in the tidal perturbation characterized by the Newtonian regres- 
sion frequency u) = 3oj 2 /4fl, which corresponds to a period of nearly 18 years 
[15,16]. It is clear from Eqs. (4)-(9) that the motion of (£} can be expressed 
as a Newtonian regression modulated by long-term (secular) relativistic pertur- 
bations characterized by the de Sitter-Fokkcr, Schiff and gravitomagnetic tidal 
effects. To illustrate this point, let us assume for the sake of simplicity that in 
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the absence of relativistic effects the lunar orbital angular momentum under- 
goes a steady regression of frequency £j and that once relativistic effects are 
included the average motion in the Fermi frame is one of precession with the 
frequency given by Eqs. (7)- (9). It then follows that the expression for (£2), 
i.e. the average of the second sidereal component of the lunar orbital angular 
momentum, contains a dominant gravitomagnetic contribution of the form 

(£ 2 )sccuiar ~ 2a/3e£(/zi?o^)wT sin(?7o + Co - w r), (10) 

where f3, ^t, Rq and Co denote, respectively, the inclination of the lunar orbit 
with respect to the ecliptic (« 5°), the mass of the Moon, the mean Earth- 
Moon separation and the longitude of the ascending node of the orbit of the 
Moon measured from the first axis of the sidereal frame. These simple con- 
siderations that are based on an initial circular orbit only indicate the nature 
of the secular terms involved; clearly, extensive calculations are necessary for a 
complete treatment. 

3 Discussion 

The results of our theoretical work are of particular interest for the description 
of dominant relativistic gravitational effects in the motion of the Moon, espe- 
cially the gravitomagnetic tidal component of the motion of the orbital angular 
momentum of the Moon. It is important to point out that the eccentricities 
of the orbit of the Earth around the Sun and the orbit of the Moon about the 
Earth should be taken into account; we have ignored them in our preliminary 
analysis [16]. As lunar laser ranging data further accumulate, it may become 
possible in the future to deduce the angular momentum of the Sun from the 
measurement of the solar gravitomagnetic contributions to the mean motions of 
the lunar node and perigee. 

It is interesting to compare our secular gravitomagnetic tidal terms with hy- 
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pothetical terms that might indicate a temporal variation of the gravitational 
"constant" G. Our results have thus far been based on a secular term propor- 
tional to r in the tidal matrix K in Eq. (3); however, as can be seen from the 
middle term in Eq. (3), similar effects could be produced if such a term appears 
in G instead. We have shown that our predictions are similar to a variation of G 
in Eq. (3) at the level of 10~ 16 yr _1 ; moreover, there are significant differences 
between the two effects that can be used to separate them [4,9]. The present up- 
per limit on \G/G\ is at the level of 10~ 12 yr" 1 ; therefore, it may be a long while 
before the gravitomagnetic effects in the motion of the Moon become detectable. 
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